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1 Introduction 



Systems with mixed classical dynamics have so far offered the most difficult problems in the 
attempts of a semiclassical quantization in terms of periodic orbits EL H, M. These problems 
have mainly three origins: 1) the existence of continuous symmetries that make (some of) 
the periodic orbits non-isolated, 2) bifurcations of stable orbits, and 3) the proximity of a 
higher symmetry that is reached by letting a continuous parameter go to zero. In all three 
cases, the original trace formula of Gutzwiller |J cannot be used because the stationary-phase 
integration transverse to the periodic orbits, that is used in its derivation, is not justified and 
leads to divergences. By now, the problems connected to 1) and 2) are essentially solved. 
Besides fully integrable systems || |J, non-integrable systems with various kinds of continuous 
symmetries can also be treated by properly extended versions of the Gutzwiller theory || 
[| . Considerable progress has also been made recently in the treatment of bifurcations, after 
earlier indications how to go beyond the simplest saddle-point integration [| [7j. Sieber and 
Schomerus [^, |S| have systematically developed uniform approximations to the most common 
types of bifurcations, expanding the action integrals in the neighborhood of a bifurcation point 



into normal forms in phase space \TU\. The resulting trace formulae interpolate continuously 
between the appropriate Gutzwiller limits that are sufficiently far away from the bifurcation 
points, where the stationary-phase integration can be applied. 

In the present paper, we shall be concerned with the third type of problems which arise 
from the breaking of a given symmetry through a continuous parameter in the Hamiltonian. 
Let us start from an integrable system described by a Hamiltonian Hq which possesses a certain 
continuous symmetry. As a consequence of this symmetry, the periodic orbits in the classical 
system occur in degenerate families living on iV-tori in phase space, where N is the number 
of degrees of freedom of the integrable system. We now perturb the system by adding to it a 
term that breaks the symmetry: 

H = H + eH 1 . (1) 

Here e is a continuous dimensionless parameter which in the following we may also call a 
"deformation" . Due to the symmetry breaking, some (or all) of the rational tori containing the 
periodic orbit families are broken up into orbits that have a lesser degree of degeneracy than 
those of H , or are completely isolated. The system ([!]) will in general exhibit mixed classical 
dynamics, and if Hi breaks all continuous symmetries of Ho, it will become chaotic for large 
values of e and for large energies, where the standard Gutzwiller trace formula can be applied. 
However, for small e the amplitudes in this formula become very large; they actually diverge 
for e — > 0. This is due to the fact that although the rational tori are broken up for e > 0, 
the periodic orbits are still not sufficiently isolated as long as their perturbed actions differ by 
less than h, and consequently the stationary-phase integration transverse to the orbits fails as 
mentioned above. One then has to find more accurate ways of performing the trace integration 
over the semiclassical Green's function; in principle, closed non-periodic orbits will thereby also 
contribute significantly to the result flTR . 

In the limit of small perturbations, e <C 1, classical perturbation theory may be used to 
derive trace formulae with finite amplitudes that yield the correct limit for e — > 0. Generalizing 
earlier attempts [O, |T3"|, Creagh fifM has recently developed a scheme to derive perturbative 



trace formulae for the breaking of arbitrary continuous symmetries, including e.g. SO (3) (spher- 
ically symmetric potentials) or SU(iV) (harmonic oscillators in TV dimensions). Applications of 
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this approach have been presented in Refs. jL5|, [16], [17], |18| . The results successfully describe 
the transitions from higher to lower (or no) symmetry for small or moderate deformations e, 
but they eventually fail when the perturbative regime is exceeded. In the limit of large pertur- 
bation, e ^> 1, one would like to recover the Gutzwiller trace formula [|lj] for isolated orbits, or 
its corresponding extension |5], |6| if some continuous symmetries are left. A closed form of an 
approximation that yields this limit as well as the correct trace formula of the integrable system 
Hq for e —>■ is called a uniform approximation, in analogy to the uniform approximations that 
interpolate continuously across bifurcations. 

Tomsovic, Grinberg, and Ullmo [[BJ have recently derived a uniform approximation for the 
breaking of U(l) symmetry in a two-dimensional system. Their result is quite general and 
applies to all systems where the rational tori are broken into pairs of stable and unstable 
isolated orbits. No analogous result is known to us for the breaking of a higher symmetry in 
any dimension. We will discuss the approach of Ref. [19| briefly in Sect. II and rederive it from 
the perturbative limit in a heuristic way that is suitable for an extension to higher symmetries 
oiH Q . 

In this paper, we derive uniform approximations to perturbed harmonic oscillators in two 
dimensions, where Hq has SU(2) symmetry (Sect. III). We furthermore apply one of the results 
to a three-dimensional cavity with small axially-symmetric quadrupole deformations, where one 
starts from SO(3) symmetry (Sect. IV). 

Our aim is not the full quantization of these systems, but the description of their gross- 
shell properties which are determined by the shortest orbits |4], [| . The use of the periodic orbit 
theory to describe shell effects in many-fermion systems in terms of a few short orbits has found 
nice applications, e.g., in nuclei, for ground-state deformations ^0] and the mass asymmetry 
of fission ||21|| ; in metal clusters, for supershells |22 and (using the perturbative approach |T4f ) 
their modifications due to deformations |IH [TH] and magnetic fields [T7| ; and in semiconductor 



quantum dots, for magnetization [13, 23 1 and conductance fluctuations [16, 24 



2 Recapitulation of U(l) breaking 



We start by recapitulating the work of Tomsovic et al. [19| for the breaking of U(l) symmetry 
in a two-dimensional system. We shall re-derive here their result in the simplified version given 
by Sieber 0, using a heuristic way which will be generalized in the later sections to systems 
with higher symmetries. 

We restrict ourselves to the most frequently occurring case that a periodic orbit family on 
a 2-torus is broken into a non-degenerate pair of stable and unstable isolated orbits. In the 
two-dimensional trace integral, one integration is performed exactly along the orbits as usual. 
The space variable q transverse to the orbit can always be mapped onto a variable <fi which is 
cyclic in [0, 2ir) and as a function of which the action shift is proportional to cos(0) ("pendulum 
mapping" in [|l!| ) . Hence the contribution of an orbit family to the trace formula is 

5 g = (e 4 * -^ r A{6)J{6) e^ SScosW d<f>\ . (2) 

I 2lT Jo J 

Here ^4(0) is the Gutzwiller amplitude function in the trace formula for the perturbed system, 
J(<p) = dq/d(j) is the Jacobian due to the variable mapping, and $ — S /h—cr 7!-/2 is the overall 
phase (including the action Sq and the Maslov index <jq) of the level density in the unperturbed 
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system with U(l) symmetry. The quantity ^S'cos (0) in the exponent is the action shift caused 
by the symmetry breaking term in the Hamiltonian. For a first inspection, perturbation theory 
might give a hint to the value of the constant 5S (which depends on the energy and on the 
parameters in the symmetry breaking term, such as deformation, nonlinearity, etc.). Hereby 
one may have to go beyond the first order of the perturbation expansion. Recent work on 



Henon-Heiles systems |1| gives us a hint that going to the lowest order at which 5S becomes 
nonzero - however high it may be - can be combined with keeping the unperturbed amplitude 
Aq. Putting A((j)) J ((f)) = Aq, which one may do in the small-perturbation limit, the <fi integral 
in (El) can be done analytically, and one obtains 

8g = A J (5S/h)cos(<S> ), (3) 

where Jq(x) is a standard cylindrical Bessel function. When 5S is zero, we have the trace 
formula in the symmetric limit (for one orbit family) 

5g = A cos($ ) = A cos(S /h - <x 7r/2) . (4) 

For large deviations from the U(l) symmetry, i.e., for 5S ^> H, we can use the asymptotic 



of Jq(x) ~ J (2/7r|x|) cos(x — tt/4) to find 



expansion 

5g ~ A ^h/27i\5S\ [cos($ + 8S/h - tt/4) + cos($ - 5S/h + tt/4)] . (5) 

This corresponds to the two isolated orbits with the action shifts ±5S and the corresponding 
corrections to the Maslov indices. Note that the two terms above arise from a saddle-point 
approximation to the integral in @ at the stationary points (fi = and </> = tt, respectively. 
The action shifts and amplitudes in (|5[) will in general be correct only in the small-perturbation 
limit. 

We want to find a uniform approximation that for 5S ^> h reaches the correct Gutzwiller 
trace formula for the pair of isolated orbits 

$9g = A u cos(S u /H - o u ii/2) + A s cos(S s /h - a s n/2) , (6) 

where the indices s and u refer to the stable and unstable orbits, respectively. We first define 
the following quantities: 

A = ~(A U + A S ), AA = ^(A U -A S ), S = ~{S U + S a ) , AS = ±(S U - S s ) , (7) 

$ = S/h-a, a = ^(a u + a s ) =a . (8) 

We now make the following ansatz for the uniform approximation, which consists of expanding 
the product ^4(0)^7(0) up to two terms with suitably chosen coefficients: 

8g u = ^2ti\AS\/Ti 3?e|e i *^^ 27r (A + AAcos0)e^ A5cosW ttyj . (9) 

For small perturbations, AS ~ SS, and therefore @ will by construction lead to the correct 
symmetric limit (^), since the divergence of the Gutzwiller amplitudes in this limit is given by 
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the first factor in Eq. (j5|). On the other hand, the coefficients in parentheses under the integral 
in (H) have been chosen such that in the asymptotic limit |Aif>| ^> H, the stationary-phase 
evaluation will lead to the amplitudes A ± AA which are precisely the Gutzwiller amplitudes 
A u and A s , respectively, and thus give the form (||). 
The integration in (|9]) can be done analytically, using 

' /- " ' <P cos <P e ix cos ^ = iJ 1 (x), (10) 
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and leads to the Tomsovic-Grinberg-Ullmo (TGU) uniform approximation [19 , in the compact 
form given by Sieber || , for the contribution of a pair of symmetry-broken isolated orbits to 
the trace formula: 

5g u = ^2n\AS\/h [A J {AS/h) cos($) - AA J^AS/h) sin($)} . (11) 

Note that this formula holds for all generic non-integrable systems in two dimensions that 
arise from an integrable system with U(l) symmetry through a symmetry-breaking term in 
the Hamiltonian that is governed by a continuous parameter. Particular examples are two- 
dimensional billiards obtained by deforming the circular billiard. The nature of the deformation 
parameter generally plays no role. The only assumption made is that the original orbit families 
(i.e., polygons in the case of the circular billiard) are broken into pairs of stable and unstable 
isolated orbits. The modification that becomes necessary when extra degeneracies due to dis- 
crete symmetries are present is trivial and will be dealt with explicitly in the examples discussed 
below. The breaking into orbit pairs is the most frequent situation. Exceptions occur, e.g., in 
billiards with octupole or hexadecapole deformations, where the boundary in polar coordinates 
is given by r(8) = R [1 + qP^cos^)] with i = 3 or 4. There the diameter orbit family breaks 
up into more than two isolated librating orbits (not counting discrete degeneracies) ||25|| . These 
have to be treated with different (and more complicated) uniform approximations. 

We also note that the deformation away from the integrable case should be small enough 
that no bifurcation of the stable isolated orbits has taken place or is about to arise. Near the 
bifurcation points, the known uniform approximations apply [0, |J which we shall not discuss 
here. 



3 Uniform approximations for SU(2) breaking 

No uniform approximation has, to our knowledge, been derived so far for systems with higher 
than U(l) symmetry. In the following, we shall do so for two systems obtained by breaking 
the SU(2) symmetry of the two-dimensional harmonic oscillator. We shall follow the heuristic 
way of deriving the uniform approximation described in the previous section, starting from the 



perturbative limit which here is treated using the approach of Creagh [14] 



For isotropic and anisotropic harmonic oscillators in any dimension, analytical trace formu- 



lae are known which converge to the exact quantum-mechanical sum of delta functions |J. [26 
For the two-dimensional isotropic case, the oscillating part of the level density is 

rS \ . 2E „ 2-kE 



r=l 
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(Note that the smooth Thomas- Fermi part is given by Aq/2.) As pointed out in [14]], the contin- 



uous degeneracy of the classical periodic orbits in this system can be described by integration 
of the surface element dQ = sin f3 d/3 d'j of a unit sphere: 

d-i r sin (3 dp = 1. (13) 
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As a result of the SU(2) symmetry, the action So is independent of the angles j3 and 7. In the 
presence of a small perturbation, the periodic orbits will be distorted, resulting in an action 
shift 5S that in general will depend on (3 and 7. Explicit ways of calculating SS(j3, 7) starting 
from a Hamiltonian of the form (]l|) are given in [|14]]. For small values of the perturbation 



parameter e, the main effects governing the level density will come from the action shift in 
the phase of the trace formula, whereas the unperturbed amplitude A can be retained. The 
perturbative trace formula for symmetry breaking then reads as 



5g pert = A Me £ M (rx) e irS °/ h \. (14) 
U=i 



Here the modulation factor M(x) (which in general is complex) is given by the average of the 
phase shift, taken over the originally degenerate periodic orbit family, 

M(x) = ^- r dj [\ml3d(3e lSS{Mln , (15) 

47T JO Jo 

and r is the repetition number. The dimensionless quantity x is proportional to e and inversely 
proportional to h, and depends on some power of the energy E. For e — * and hence x — > 0, 
we have SS — > so that M —> 1, and the unperturbed trace formula ([12]) is recovered. The 
repetition number r in (|14]) cannot be summed up to arbitrarily high values, since the argument 
rx must remain of order unity or smaller for the perturbation approach to be valid. Hence, the 
maximum value r m must be chosen such that r m x < 1 for given values of e and E. This excludes 
in general the possibility of quantizing the system through the trace formula in this approach. 
However, we shall be interested only in the low-frequency components of the oscillating level 
density, i.e., in its gross-shell structure that is governed by the shortest periodic orbits and their 
first few repetitions [J. 

Our main task now is to generalize the modulation factor (]T5|) in such a way that the trace 
formula ( |T4] ) goes over to the Gutzwiller formula Jl| for isolated orbits in the limit of large 
perturbations e that fully break the symmetry, whereas the limit yVl(O) = 1 is preserved. If we 
succeed in finding such a generalization, it will smoothly interpolate between the exact trace 
formula (0) for the harmonic oscillator and the Gutzwiller formula for the symmetry-broken 
limit, and hence be a suitable uniform approximation. Note that Eq. ( |T4"P with flT5| ) is exactly 



of the same form as Eq. (|2|) for the U(l) case, except that we now have a two-fold integral. We 
can therefore take the same course of action to find a uniform approximation that has the above 
two limits: (i) evaluate the asymptotic amplitudes for large values of x (i.e., for large e), (ii) 
map the exact action shift 5S onto the form obtained from perturbation considerations (but 
with freely adjustable parameters), (Hi) include under the integral a parameterized expression 
of the same form for the product of the amplitude function A and the Jacobian J of the 
mapping function, and (iv) adjust all the parameters such that in the asymptotic limit of large 
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x the correct Gutzwiller amplitudes and actions of the isolated orbits are obtained, while the 
unperturbed limit fll2|) is kept for x — 0. 

This procedure is very similar in spirit to that used by Sieber and Schomerus in their uniform 
treatment of many types of bifurcations || . However, instead of starting from a phase-space 
representation of the trace formula and expanding the action in normal forms, we use the group 
integral representation of the (unperturbed) trace formula developed by Creagh and Littlejohn 
0. The latter exploits directly the properties of the symmetry group characteristic of the 
system Hq] it has been used as a starting point for the perturbative trace formula of Creagh 
rjJJ, as shown above in Eqs. (JT^fT^). 

Different from the U(l) case (with the exceptions mentioned at the end of the previous 
section), the explicit results which we obtain here do depend on the explicit form of the Hamil- 
tonian Hi. On the other hand, one of the results (that for the standard Henon-Heiles potential) 
turns out to apply also to a three-dimensional cavity with small axially-symmetric quadrupole 
deformations. This is due to the close relation between the SU(2) and the SO(3) symmetry 
that is broken in the latter case. We will discuss this system in section IV. 

3.1 The quartic Henon-Heiles potential 

We first investigate the quartic Henon-Heiles (HH4) potential in two dimensions, given in polar 
coordinates by 

V(r,9) = icuV-|r 4 cos(40), (16) 

which has recently been investigated [0 in the framework of the semiclassical perturbation 
theory described above. The limit a = is a harmonic oscillator with the SU(2) symmetry. The 
anharmonic term makes the system non-integrable with mixed classical dynamics. It retains 
a discrete four-fold rotational symmetry with four saddle-points at the energy E* = uj A /Aa, 
through which the particle can escape. The shortest periodic orbits in this system are two pairs 
of straight-line librating orbits (named Ai and A2 in Ref. [|H|) and a circulating orbit (named 
C). Like in the standard cubic Henon-Heiles (HH) potential |18|, [27], 28 1, the system is scaled 



with a and its dynamics can be described in terms of one continuous parameter, the scaled 
energy e defined by 

e = E/E* = AaE/uj\ (17) 
The actions of the periodic orbits are changed in first order of the perturbation by a shift 

3 

SSi = hx (sin 2 /? cos 2 7 — cos 2 /?) , Tlx = S — e . (18) 

Then, the perturbative modulation factor J\4(x) in fll5[) becomes (with u = cos/3) 

M(x) = ^- d-f J\m(3d/3e lx{sin ^ cos ^~ cos2 ^ = ^ j*™ d 1 J* du e *=[(i-« 2 )c°s 2 7 -u 2 ]_ (19) 

An alternative expression of this is obtained by rotating the unit sphere about an angle 7r/4 
along the 2 axis: 

11/ ,11 

ei = -j= ei - -j= e 3 , e 2 = e 2 , e 3 = -j= e a + -j= e 3 , (20) 
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which leads to 



1 



I" sin (3 dp e ix sin(2/3) COS7 . 
Jo 

This integral can be performed analytically (see |18|]) and yields 



Jl/4 



X 



(21) 



(22) 



The asymptotic expansion of the Bessel functions J M (x) for large x leads to three terms which 
we label according to the names of the periodic orbits: 



M(x) ~ M c {x) + M Al (x) + Ma 2 {x) 



(x > 1) 



(23) 



with 



M c (x) 
M Al (x) 
M A2 (x) 



1 

2x 



2y/2x 
1 

2^2x 



cos(+x — 7r/2) 
cos(— x + 7r/2) 



A c ~ — 



Ao 



2y/2x' 
_j4o_ 

2^0;' 



(24) 



The first term corresponds to the loop orbit C, the second and third terms to straight-line 
libration orbits Ai and A2, respectively; each of them has a discrete degeneracy of two which 
is included in the amplitudes. These are the only periodic orbits with periods of order 
To = 2ir/u up to energies e ^ 0.85. The forms on the left-hand side of Eq. (|24]) contain the 
exact Maslov indices ac = 0, oa x = +1, and oa 2 = — 1 of the isolated orbits and, with Eq. 
([12"|), yield the asymptotic amplitudes of the trace formula shown on the right-hand side of Eq. 
(|24|). For small x, these amplitudes have been shown in Ref. [18| to reproduce numerically 
well the diverging Gutzwiller amplitudes of the isolated periodic orbits. Note that the action 
shifts predicted at the first order of the perturbation theory are zero for the orbit C and ±hx 
for the orbits Ai and A2, respectively. This can be checked numerically for all orbits, as well 
as analytically for the orbits Ai and A2. The actions of the latter, being straight-line one- 
dimensional integrals, can be expressed analytically in terms of complete elliptic integrals (see 
the appendix) and then be Taylor expanded in powers of e. The result is 



Sa 1 



Sa 2 



So(l + -e + 
Sod-|e + 



35 

1024 

35 

1024 f 



e 2 + 



1155 

65536 
1155 

65536 



e J + . . . 



e A + . . . 



(25) 



which confirms the first-order action shifts ±hx given in Eqs. ( |r8| , |2"4j ) . 

For the following, it is important to trace back the asymptotic forms ( p4|) to singular points 
of the integral (|19"|). For this, we first evaluate the asymptotic form of the u integral. The 
stationary point at u = yields a term ~ l/^/x, and the end point at u — 1 yields a term 
~ 1/x. The latter, after an exact integration over 7, gives the contribution AiA 2 {x)- The 7 
integral over the first term has four stationary points; saddle-point integration at 7 = and 7r 
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yields the contribution J\4a 1 (x), and at 7 = n/2 and 37r/2 it yields the contribution J\A.c{x). 
In summary, we get the periodic orbit contributions asymptotically as follows: 



orbit A 2 
orbit Ai 
orbit C 



from u = l, 
from u — 0, 
from u = 0, 



(any 7) , 

7 = and it , 

7 = 7r/2 and 3tc/2 



(26) 



We now construct a uniform approximation which for large x yields the correct asymptotic 
Gutzwiller trace formula with the contributions from the three leading isolated orbits A 1; A 2 , 
and C. For describing the gross-shell structure of the level density, it is sufficient to include 
only the lowest few harmonics,]] i.e., the first few repetitions r of the primitive orbits, in the 
trace formula. In the following, we shall only give the results for r = 1; higher repetitions can 
be included according to Eq. fl!4D . Hence, the Gutzwiller limit is written as 



8g, 



G 



A Ax cos 



S Ax n \ , 

— -^2 ) + a a 2 cos 



S 



A 2 



h 



n \ , A ( Sc 71 

aA *2) +AcC ° S V%- ac 2 



(27) 



This limit can be imposed by including under the integral fljjp the product of amplitude func- 
tion A(u, 7) times Jacobian ^(tt, 7) in the same form as in the exponent, but with different 
parameters: 



M u {x) = — [ % dj f d ue ^[(i-" 2 )<*» 2 7-u 2 ] 
2tt Jo Jo 



2x — — 

x— [\p2A A2 u 2 + (1 - u 2 ) [V2A Al cos 2 7 + A c sm 2 ^\ } .(28 



A 

Note that the coefficients in the second line are precisely the inverses of the asymptotic ampli- 
tudes given in In this way, the approximation ( |2"8"D leads by construction to the Gutzwiller 
limit (|27D for large x, whereas for x — > the diverging Gutzwiller amplitudes exactly cancel 
altogether and M. u — > 1, as required. 

The integrals occurring in (p8|) can all be done analytically. For the integral appearing as 
the coefficient of Ac, we note that e 2 , = sin 2 /?sin 2 7 is invariant under the rotation (|20|) , so 
that we can replace the phase in the exponent by that appearing in Eq. (pl|). Then, using the 
same transformations as in Ref. [18| for obtaining the form (^) of the perturbative modulation 
factor, and exploiting some recurrence relations amongst the Bessel functions J M (x), we obtain 

— f 277 d 7 - ^ sin p dp sin 2 /? sin 2 7 e ix sin{2/3) COS7 
2ir Jo 2 Jo 



71 



4V2 



J -V4 ( o ) J l/4 



3/4 



X 



J: 



3/4 



(29) 



The other two integrals can be found first by taking the derivative of Eq. fll9|) with respect to 
x, and second by integrating the terms proportional to u 2 in (|28|) over m by parts, and finally by 



*A more consistent truncation of the trace formula, followed below, is achieved by Gaussian smoothing over 
a small energy range 7, resulting in the amplitudes multiplied by exponential damping factors; see Ref. |p8| for 
details. These exponential factors are included in the amplitudes A Al , Aa 2 , and Ac and their combinations 
defined below. 
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taking suitable linear combinations of the results of these two operations. The final expression 
is 

9/r r _ _ _ 

(30) 



M u {x) = — A c M-{x) + V2A A M+(x) - iV2~AA A M'{x) 



where 



and 



A 



A 



A Al + A M 



AA. 



A, 



A M 



(31) 



M±(x) 
M\x) 



71 



4V2 

7T 



J_ 1/4 ( f ) Jl/4 (f 



± J_ 



3/4 



3/4 



4^ 



J_ 



1/4 



I) J5/4 (I) + ^ (I) J3/4 (I 



(32) 



The last step now is to insert the modulation factor ( |30"D into Eq. (pf ) and to replace the 
unperturbed action So in the phase of the first term of (|30|) by Scs and for the other two terms 
by the average action S A of the orbits Sa 1 and S A2 , while the perturbative action shift hx is 
redefined as their difference A Sa 



hx = AS A = - (S Al - Sa 2 



$a = - {Sa 1 + Sa 2 



(33) 



By these replacements, we ensure that the phases in the asymptotic level density (|27D contain 
the correct numerical actions of the isolated orbits. 

The final form of the uniform level density for the HH4 potential (including only the prim- 
itive orbits, i.e., r = 1) is then: 

5g u = (2AS A /h)^e{e^oA c M-(AS A /h) 

+V2e% 3A [A A M+ (AS A /h) - % AA A M' (AS A /h)] } 
= (2AS A /h) {A C M- (AS A /h) cos(S c /h) 

+V2[A A M + (AS A /h) cos(S A /h) + AA A M' (AS A /h) sm(s A /h)]}.(34) 

As in the case of the TGU uniform approximation discussed in Sec. II, this formula only 
holds as long as no bifurcations of the stable isolated orbits occur. For the primitive orbits of the 
HH4 system, this does not happen up to e ~ 0.85, where the primitive orbit Ai undergoes an 
isochronous bifurcation. In the numerical results shown below, we have smoothed the Gutzwiller 



amplitude of this orbit (see Ref. [18] for details) in order to simulate a better treatment of this 
bifurcation by the corresponding uniform approximation ||. 

In Fig. we show a numerical compilation of the three semiclassical approximations to the 
level density 5g discussed here (given by dashed lines) and compare them to the quantum- 
mechanical one (given by the solid lines). In all cases, Gaussian averaging over the energy 
E with a range 7 = 0.6 hu has been applied. This damps the amplitudes strongly enough 
so that only the primitive orbits (r = 1) need be kept in the semiclassical approaches. At 
the bottom of this figure, the Gutzwiller result (|2~7D for the isolated orbits is shown. It gives 
an excellent agreement with the quantum result for E > 10 hu, even up to e ~ 1 (the saddle 
energy is E* ~ 39 hu for the case a = 0.0064 chosen here). For small energies, it diverges 
due to the approaching SU(2) symmetry limit. In the middle of the figure, the perturbative 



10 



result is shown. It reproduces the quantum result up to E ~ 13 hu, thus catching the 
essential feature of the symmetry breaking. (Note that at E ~ 13 %u we have x ~ 3, showing 
that the perturbative approach may be used also for values of x somewhat larger than unity, 
cf. [[L4|, |i~8fl.) For E > 14 hu the anharmonicity is, however too large for perturbation theory to 
apply. Finally, the top part of Fig. [I] shows the present uniform approximation that reproduces 
the quantum result at all energies. 
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Figure 1: Oscillating part of level density of the HH4 potential (a = 0.0064), averaged 
over a range 7 = 0.6, versus energy E (units: huj). Solid lines: quantum-mechanical 
results. Dashed lines: semiclassical approximations (with r m = 1). Bottom: Gutzwiller 
trace formula (diverging at small energies); middle: perturbative trace formula by Creagh 
(failing at large energies); top: present uniform approximation (working at all energies). 



In Fig. I we compare the results of the present uniform approximation (dashed lines) with 
the quantum results (solid lines) for two different values of the anharmonicity parameter a. In 
the two upper panels, the same smoothing width 7 = 0.6 %uj was used as in the previous figure 
and only the primitive orbits (r = 1) Ai, A2 and C were included. In the two lower panels, 
we have used a smaller smoothing width 7 = 0.25 fro; and included the second repetitions (i.e., 
r m = 2) of the three orbits as well. Now a more detailed fine structure of the level spectrum is 
resolved; even this is well reproduced by the semiclassical approximation. Only in the regions 
corresponding to e = E/E* > 0.8 (E ^ 31 hu), some small differences can be seen which are 
explained by the fact that at these energies, new orbits with actions comparable to those of 
the included orbits with r = 2 exist (after period-doubling bifurcations of the orbits C and Ai) 
which we have not included. Their inclusion would necessitate a proper uniform treatment of 
the corresponding bifurcations, which is outside the scope of this paper. 

In the actual computation of the uniform trace formula for small energies e, cancellations 
between the diverging Gutzwiller amplitudes take place. This requires their rather accurate 
numerical determination. For practical purposes, it is advisable to take advantage of the fact 
that for sufficiently small arguments x = ASa/^, the result ( ^4[) goes over into the perturbative 
trace formula ( JH| ) with the modulation factor ( |2"2| ) which is numerically much more robust. We 
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found that one may switch between the two formulae for values 1.5 ^ x ^ 2.5 without visibly 
changing the results shown above. 
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Figure 2: Oscillating part of level density of the HH4 potential versus energy E (units: Too). 
Solid lines: quantum-mechanical results. Dashed lines: semiclassical results with present uniform 
approximation. Left side: a = 0.0036 (E* ~ 69?k<j), right side: a = 0.0064 (E* ~ 39hui). Upper 
part: 7 = 0.6 and r m = 1, lower part: 7 = 0.25 and r m = 2. 



3.2 The standard Henon-Heiles potential 



We will now investigate the original Henon-Heiles (HH) potential |27 , given in polar coordinates 
by 

1 rv 

V(r, 9) = - 00 V - - r 3 cos(30) . (35) 

2 3 



The three shortest isolated orbits in this potential, which have been shown [f28| to govern the 
beating gross-shell structure in the level density of this system, are a straight-line librating orbit 
A, a curved librating orbit B, and a circulating orbit C (similar to that in the HH4 potential). 
Orbits A and B have a discrete degeneracy of three, due to the three-fold discrete rotational 
symmetry of the potential (^), and the orbit C is two- fold degenerate due to time reversal 
symmetry. 

Because of the odd power of the anharmonic term in ([35]), the first-order perturbation result 
is zero both classically and quantum-mechanically. In second-order perturbation, the action 



shift SS due to the anharmonicity is given [18] by 



5S 2 = frx(5-7cos 2 /?)/6, (36) 
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with 

S E S u 

fix = = — e, E = — -. 37 

12 E* 12 6a 2 ( J 

Here E* is the saddle energy for the potential (|35|) , and So is the action of the unperturbed 
harmonic oscillator given in (|i~2"l). The perturbative modulation factor A4(x) becomes 

M(x) = ^JdSl e »( 5 - 7 " 2 )/ 6 = £ du e -( 5 - 7 " 2 )/ 6 . (38) 

This integral can be expressed analytically in terms of the Fresnel functions C(x) and S(x) (we 
use the convention of Abramowitz and Stegun [29]): 



M(x) = e 5 -/ 6 1 [C(0 - iS(0] , e = y ^ ■ (39) 
Using the asymptotic expansion of the Fresnel functions for large arguments x ^> 1, 

C(0 ~ ^ + 4 sin KV2), 
2 7r4 

5(0 ~ 5 - ± cosK 2 /2) , (40) 

we get the asymptotic result 



/ 3 3 
^ ~ \ 77Tl A v c°<So/h + 5x/Q - tt/4) + cos(5 //i - x/3 + tt/2) , (41) 

y 14|x| 7\x\ 

which is again only correct for sufficiently small values of the scaled energy e so that x 1. 
Note that (|4lD predicts the action shift ASc and the average action shift ASab to be 



AS C = S C -S = -hx/3 , ASab = S AB - S = 5hx/6 , Sab = \(Sa + S b ) • (42) 



This is numerically well fulfilled at lower energies. The expansion (|4T| ) is also obtained from the 
asymptotic analysis of the integral fl38D , whereby the first term comes from the stationary point 
at u = and the second term is the end-point correction from u = 1. The latter corresponds 
to the C orbit in the HH potential with the correct Maslov index [18 ac — — 1- Its amplitude 
is thus predicted in the second-order perturbation theory as 

pa 6E 36 I 

c 7(hu) 2 x 7nhu e 1 1 



which was shown in Ref. |L8| to describe well the numerical Gutzwiller amplitude Ac, not only 
for small e where it diverges, but also up to energies e ~ 1. The first term in (41) corresponds 
to the sum of orbits A and B which, however, are still degenerate on the circle 7 e [0, 2tt) at 
the second order of the perturbation expansion and therefore have an asymptotic amplitude 
too large by a factor 1/y/h. These orbits are broken up only at fourth order, with an action 
shift 

5S 4 = -%sin 3 /?cos(37) , y = c 4 e 3 . (44) 
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Including this into the phase of the modulation factor, we get 

M(x,y) = / due ix ^- 7u2)/6 — d 7e -«/(i- 2 ) 3/2 cos(3 7 ) 
Jo 2ir Jo 



o 



due^- 7u ^ 6 J [y(l-u 2 f /2 ]- (45) 



In the asymptotic expansion of this modulation factor, the Bessel function modifies only the 
contribution from the stationary point at u = 0, leading to 



/ Q Q 

5 9 ~ \ TJn A Mv) ^s(S /h + 5x/Q - tt/4) + —A C os(S /% - x/3 + tt/2) . (46) 
y 14|x| 7|x| 

Upon further expansion of Jo(y) for y ^> 1, this yields the amplitudes and the correct Maslov 
phases [[TH] = 1 and erg = of the (now isolated) orbits A and B. The amplitudes, which are 
equal at this order of the perturbation theory, go like 1 / \J\xy\ and thus have the same power 

of % as Aq 2 . The actions of the orbits A and B are shifted from the average value Sab by an 
amount ±hy, respectively. For the orbit A, this can again be checked by analytical integration 
(see the appendix) and Taylor expansion of its period. The action becomes 

_ /, 5 385 2 85085 3 \ 
^ = M 1+ 72 e+ 15552 e + 6718464 e + ' ' ') " (4?) 

The first correction term is the average action shift ASab (0)? obtained for both orbits A and 
B at second order of the perturbation theory; the next term gives the value of C4 in Eq. 
With this, the fourth-order prediction of the average amplitude of orbits A and B becomes 



- E J 3 216 J 3 J- U8) 
AB (hu)^7\xy\ 77r^V55e 3 / 2, { } 

[Note that this amplitude contains the degeneracy factor 3 of the orbits A and B mentioned 
above; it is related to the factor 3 in the argument of the cos function in (44). Similarly, A^j 2 
(43) contains the time-reversal factor 2 of the orbit C, which is related to the two end-point 
corrections, coming from the u integral in (|38|) which originally runs from —1 to +1, each giving 
one half of the second term in (Wm!-\ 

We now redefine the quantities x and y in terms of the true actions of the isolated periodic 
orbits 

hy = 5S = ^(S A - S B ) , ^x = AS = Us A + S B ) - S c , (49) 
I b 2 

and introduce the combinations of the Gutzwiller amplitudes 

A AB = \ {A a + A B ) , AA AB = g (Aa ~ A B ) . (50) 
The uniform approximation to the modulation factor is then 

AqJo 2-kJo 

2\AS\ . 2 [a\AS\ 2 , l27c\5S\ v- ... 

— Ac m + y (1 — M )y — ^ — \A-ab + AAab cos(37) }. (51) 
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The second line above is again the parameterized product A(u, 7) J^(u, 7), with coefficients 
chosen such that asymptotically for large AS" and SS, the Gutzwiller limit 



5gc = Aa cos(Sa/Ti — &ATr/2) + A B cos(Sb/Ti — o- B n/2) + A c cos(Sc/h — acn/2) 



(52) 



is reached. (The degeneracy factors due to discrete symmetries discussed above are again 
included in the amplitudes.) In the limit AS* = 5S = 0, on the other hand, M. u (|5T| ) still 
reduces to unity as it should. 

The two-dimensional integrals in (^T|) cannot be done analytically here. It turns out, how- 
ever, that we can approximate them without violating the above two limits.Q For that we 
note that the first term (corresponding to orbit C) asymptotically gets contributions only from 
u ~ 1; we may therefore put u — 1 in the exponent of the 7 integral which then becomes 
unity. The second term of (|5lD , corresponding to the orbits A and B, asymptotically only 
gets contributions from u ~ 0. Putting u = in the exponents of the 7 integrals leads to 
Bessel functions Jo an d Ji like in the TGU formula (|H|) . The remaining u integrals can now 
be expressed again in terms of the Fresnel functions after some partial integrations; hereby we 
keep only the leading-order terms in h. We then arrive at the uniform trace formula for the 
HH potential, including the contributions from the three primitive orbits A, B, and C, 

Sgu = A c cos(- + -) 



+ C 



+ 



An\5S\ 




(S c ^ 



A , T (5S\ . [S c 71 



(53) 



In the low-energy limit AS — > 0, 55' — ^ 0, the second line of (|53|) cancels the orbit C term in 
the first line, the fourth line vanishes, and the third line yields the HO trace formula 
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At energies where the C orbit is well isolated and separated from the A and B orbits 
(AS ^> h) but the splitting of A and B is still small (SS < h), we can use the asymptotic forms 
POD of the Fresnel functions. Hereby the second terms from (fi0| ) combine to cancel the second 
line in (|53D, whereas the leading terms combine into 



Sg u 



.' Sc 7T , 

A c cos( T + -)H 



Aab Jo I -z- I cos 




(x > 1J54) 



so that the second line contains nothing but the TGU approach to orbits A and B, kept separate 
from the isolated orbit C contribution. 



^These approximations correspond, in fact, to neglecting only higher-order terms in Ti. We have checked by 
numerical integration that these do not affect the results discussed below. 
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In Fig. H we show the numerical results obtained with the uniform approximation (|53|) using 
r m = 2 (dashed line), compared with the quantum-mechanical result (solid line) for the HH 
potential with a = 0.04 and a Gaussian averaging with width 7 = 0.25 hu. The agreement 
is perfect for all energies up to e ~ 0.73. The differences seen for e > 0.73 are due to some 
missing orbits which arise after period-doubling bifurcations of the primitive orbits A and C, 
and perhaps to the isochronous bifurcation of the A orbit that occurs only at e ~ 0.97 but 
may make itself felt in the amplitude A a already at smaller energies. (Note that, in contrast 
to the HH4 case above, we have not smoothed this amplitude here.) All these bifurcations can 
be treated with the uniform approximations developed in Refs. ||. Part of the disagreement 
for e > 0.75 is also due to inaccuracies in the diagonalisation of the quantum HH Hamiltonian 
in a finite basis 
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Like in the case of the HH4 potential, the SU(2) limit is reached in the uniform approxima- 
tion (|53"D through a subtle cancellation of divergences in the Gutzwiller amplitudes. Numerically, 
the most robust procedure is to use the perturbative result (^) for values of x — 6AS/7H up 
to ~ 1.5 — 2.5, and then to switch to the uniform approximation. 





0.45 0.5 0.55 0.6 0.65 0.7 0.75 

Figure 3: Oscillating part of level density of the Henon-Heiles potential (a = 0.04), Gaussian 
averaged over a range 7 = 0.25, versus scaled energy e. Solid line: quantum-mechanical result; 
dashed line: semiclassical result in the present uniform approximation with r m = 2. (Units: 

flLJ.) 
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4 Axial quadrupole deformations of a spherical cavity 



We finally discuss three-dimensional cavities with axially symmetric quadrupole deformations 
that can be treated with the uniform approximation developed above for the HH potential. 
The boundary of the cavities is in polar coordinates given by 

R(9,(P) = Ro[l + eP 2 (cos9)], (55) 

where P2{x) = (3x 2 — 1)/2 is the second Legendre polynomial. The periodic orbits in a spherical 
cavity with ideally reflecting walls have been discussed extensively by Balian and Bloch , and 
an analytical trace formula has been given by these authors for its level density which we again 
write as in Eq. @. The three-dimensional degeneracy of the polygonal orbits with N > 3 
corners, due to the SO (3) symmetry of the sphere, can be described by the three Euler angles 



(a,/3,7), and the corresponding group integral is || [14 

— 7T / dq I sin (3 d/3 da = 1 . (56) 
8n z Jo Jo Jo 

We choose the Euler angles such that f3 is the angle between the normal to the plane of motion 
of an orbit family and the z axis, 7 describes the orientation of a single orbit within this plane 
(i.e., within the family), and a describes rotations of the orbit plane around the z axis. As 
long as we restrict the deformations of the cavity to be axially symmetric (around the z axis), 
the action of the perturbed orbits will not depend on the angle a. For small deformations, 
the contribution of each periodic orbit family to the spherical trace formula can thus again be 
written in the perturbative approach as 

5g pert = {M{x, y ) e ^ n -^} , (57) 

where the modulation factor is 

1 i-2-k r l 

M(x,y) = — dj due i5s{u ^ y) . (58) 
2rr Jo Jo 

Here x and y are parameters depending on the deformation e, the energy, and the specific orbit 
type. When both x and y are non-zero, the orbit families are broken into singly degenerate 
families of orbits lying in planes. These come in two types. One type is lying in planes 
perpendicular to the symmetry axis (coming from the end-point correction at u — 1); we will 
henceforth call them the equatorial orbits. They are identical to the orbits in a two-dimensional 
circular billiard. The other type of orbits lie in planes that contain the symmetry axis (coming 
from the stationary point u = 0); we will call these the planar orbits. They are isolated in their 
plane of motion and will, for small deformations, come in pairs of stable and unstable orbits. 
With the exception of a completely isolated linear orbit that oscillates along the symmetry 
axis, all the equatorial and planar orbits have the continuous one-parameter degeneracy due to 
rotation of the plane of motion around the symmetry axis. 

We see therefore that the situation is completely analogous to that of the HH potential 
discussed in the previous section: orbit C corresponds to the equatorial orbit families (with two 
orientations connected to time reversal), and the orbits A and B to the pairs of planar orbits. 
The only difference is that, since we start here from a three-dimensional system with SO(3) 
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symmetry, the breaking of two degrees of degeneracy leaves the orbits in singly degenerate fam- 
ilies. Furthermore, the planar orbits do not have the discrete three-fold rotational degeneracy 
as the A and B orbits in the HH potential, but occasional two-fold degeneracies due to the 
symmetry of reflection at the equator plane. We can therefore take over the previous results 
with little effort. We shall first discuss the family of polygon orbits with N > 3 corners, and 
come to the diameter orbits (N = 2) later. 



4.1 Polygon orbits with N > 3 

The families of polygon orbits in the sphere have in principle the full SO (3) symmetry. However, 
when we limit the deformation of the cavity to axially- symmetric shapes, the action shift 
depends only on two of the Euler angles, so that the symmetry breaking problem is the same as 
in the HH potential considered above. In first-order perturbation theory [pL5| | , the action shift 
depends only on the angle (3: 

8 Si = -AS P 2 (cos (3), AS = ^S . (59) 

This corresponds to the fact that a quadrupole deformation to lowest order in e is identical 
to that of an axially symmetric ellipsoid, which is an integrable system and thus has an extra 
symmetry with two-fold degenerate orbit families (i.e., polygons with N > 3 fitting into the 
ellipse that gives the boundary of the plane of motion). As for the orbits A and B of the HH 
system, this symmetry will be broken only at a higher order in the perturbation expansion; in 
the present case we expect this to happen at the third order in e, with an action shift 

5 S3 = —5S p i sin 3 (3 cos(n7) , (60) 

where n counts the discrete symmetries of the planar orbits and will again be absorbed into 
the Gutzwiller amplitudes. With this, the modulation factor becomes 

M (x, y) = f 1 du e^ 1 " 3 " 2 )/ 2 — f^d-y e ~^-^? 12 °°»(m0 ( 6 l) 



lo 2tt Jo 

with Tix = AS and hy = 8S p i. For small x and y, the action shifts will be 

(u = 1) : S eq — So = —x = —HAS , 

(u = 0) : S p i -S = +x/2 = +KAS/2 , (62) 

with 

S p i = —(Spi + Spi) , Sp t = Spi + 8S p i j = Spi — 5 Spi . (63) 

The modulation factor ( |BT] ) has exactly the same form as that of the HH case (|£5|). However, 
different from there, x and y can have both signs, depending on the sign of e. Therefore, 
switching from prolate deformations (e > 0), which is analogous to the HH situation, to oblate 
deformations (e < 0), one has to take the complex conjugate of the modulation factor before 
using it in (|57|) . This only affects some of the signs in the final results; in the following formulae 
the upper signs correspond to the prolate case and the lower signs to the oblate case. 
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The final trace formula for one type of orbits (with fixed number N > 3 of corners, and 
without including their higher repetitions) becomes 



(JV>3) 



A eq COS $ e(? 



CI 



± s\ 



\2\AS\ 
f 2\AS\ 



2\5S\ 
\]AS~\ 



I4tt\5S\ 



IAtt\5S\ 

h 



(6S\ 
— J COS $ eg - 

Apl Jo I I COS $ p ; 



. , T (5S\ . ^ 
AA pZ Ji I — I sin $ e? 

- AA pl J x ( — \ sin $ pi 



A p i J 1 -r- I sm $ pZ + Aipj <A ( -H cos ®p 1 



(64) 



where we have redefined 



— Tlx — AS* — Spi S e q , 



h V = 8S=\{S£-S$ t 



5, 



7T 



7T 



(7 - 



and the combinations of Gutzwiller amplitudes (including discrete degeneracy factors) 



A p i — -^(Apf + A s pl ) , AApi — -^(.Api — A 
The correct Maslov indices of the asymptotic orbits become 



eq 



ao =F 1 , Opi = ]- (a pl + a s pl ) = a ± -, a pl = a s pl + 1. 



(65) 



(66) 



(67) 



At moderately large deformations, the equatorial and planar orbits are sufficiently well 
separated, so that we can take the limit \AS\ Ti. The result (pi) then simplifies to 



(N>3) 



A eq cos $ e(? 



+ 



/ 2tt|(5S'| 



(68) 



which again contains the simple TGU formula applied to the planar orbits. At larger defor- 
mations (but still before any bifurcations take place), we expand the above for I^S*! ^> Ti and 
obtain the Gutzwiller limit 



A eq cos $ eo + Al v cos + A% cos 



eq 



pi 



-pi- 



rn 



4.2 Diameter orbits (N = 2) 

The diameter orbits (N = 2) in a spherical cavity have only a continuous degeneracy of two, 
since rotation about themselves does not bring about any new orbit 0. It is convenient to 
redefine the Euler angles such that (3 describes the angle between the diameter orbit and the 
z axis. In first-order perturbation theory, the action shift of the diameter orbit due to a 
quadrupole deformation (^) is then given by [T5|] 



SSi = eSo-P2(cos f3) 



(70) 
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Hence the modulation factor for small deformations becomes 



Mix) 



due ix^-l)/2 



(71) 



with Tix = eSo. Asymptotically, we get from u = the equator orbits which keep the U(l) 
degeneracy corresponding to rotation about the z axis, and from the end point u = 1 we find 
the isolated diameter orbit along the z axis. Their action shifts for small x will be 



Seq — So — —fix/2 , Si so — Sq — +?IX 



(72) 



The situation thus corresponds exactly to the case for the polygonal orbits, but with y = since 
there is no further splitting of the isolated diameter orbit, and with the roles of prolate and 
oblate deformations interchanged. The uniform trace formula for the diameter contributions is 
thus 



(JV=2) 



+ 



V2A 



(:'</ 



CI 



'2\AS\ 



COS 



S 




Hereby we have (re) defined 



AS = S lso - S eg , $ 



eg 

n 



Siso | -, 

—Z CisoTT, O'iso — C"0 ± - 1 - 

n 2 



7T 



^0" 



.(73) 



For \AS\ ^> h, we get the Gutzwiller limit for the diameter orbits: 



A eq cos $ eg + A iso COS $ is 



with 



S, 



71 



eq 



- a, 



eq 



' eq 



o"o =F 



(74) 



(75) 



(76) 



In all equations above, the upper and lower signs correspond to AS > (prolate case) and 
AS < (oblate case), respectively. 

In Fig. |] we show the oscillating part 5g(k) of the level density for a very small quadrupole 
deformation of e = 0.01, Gaussian-averaged over k with a range 7 = 0.6 Rr 1 and plotted versus 
the wave number k = \ / 2mE/h. Note the pronounced "supershell" oscillations that are mainly 
a result of the interfering triangle and square orbits 0. In the semiclassical result (dashed 
line), periodic orbits with up to iV = 6 reflections were included in the uniform approximation 
given by Eqs. fl6~4],|73D ; identical results are obtained at this deformation in the perturbative 
approach. Both reproduce very well the quantum- mechanical result (solid line). For details 
concerning the calculation of the Gutzwiller amplitudes for the equatorial and planar orbits in 
non-integrable axially deformed cavities, we refer to a forthcoming publication 
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Figure 4: Oscillating part of level density in an axially-symmetric quadrupole-deformed 
cavity with deformation e = 0.01, Gaussian averaged over k with a range 7 = 0.6, versus 
wave length k (units: R^ 1 ). Solid line: quantum-mechanical result; dashed line: semiclas- 
sical result in the present uniform approximation using equatorial and planar orbits with 
up to N = 6 reflections. 



In Fig. H], we show the corresponding results at a quadrupole deformation of e = 0.1 at 



which the supershell beating is already reduced [HJ] 



quantum 
uniform 
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kR 

Figure 5: Same as Fig. 0, but for the deformation e = 0.1. 
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The relative importance of the different orbits at these deformations can most easily be 
analyzed in the Fourier spectra of the oscillating level density 5g(k). Since billiard systems 
exhibit scaling (i.e., the properties of the periodic orbits do not depend on energy), the Fourier 
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L/R 

Figure 6: Fourier transform of the level densities shown in Fig. |4j. The symbols indicate the 
periodic orbits corresponding to the Fourier peaks. The deformation (e = 0.01) is so small 
here that equatorial and planar orbits are not separated yet. The perturbation approach 
(short-dashed line) gives practically identical results as the uniform approximation (long- 
dashed line) and the quantum mechanics (solid line). 




3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 

L/R 

Figure 7: Fourier transform of the level densities shown in Fig. 0. At this deformation 
(e = 0.1), the equatorial orbits (symbols in heavy circles) are well separated from the planar 
orbits (symbols in thin ovals). However, the pairs of isolated stable and unstable planar 
orbits are not separated yet, as seen from the large peaks obtained with the Gutzwiller 
trace formula (thin-dashed line). The tiny bump at L ~ 4.4 R corresponds to the unstable 
isolated diameter orbit along the symmetry axis, whose Gutzwiller amplitude is smaller 
than all the others by a factor y/KJkR. 

transform with respect to k gives directly power spectra in which the peaks occur at the 
lengths of the contributing periodic orbits. Figures ^] and |7| present the absolute values of 
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the Fourier transforms of the level densities shown in the above two figures. In Fig. || the 
peaks corresponding to equatorial and planar orbits cannot be separated, since at the small 
deformation e = 0.01 the spherical tori are hardly broken. Indeed, the perturbative result 
(short-dashed line) gives here practically identical results as the uniform approximation (long- 
dashed line); both agree very well with the quantum result. At the deformation e = 0.1, the 
equatorial and planar orbits are already well separated, as can be seen from Fig. [7], so that the 
full uniform approximation flS3|) for the orbits with N > 3 can here be replaced by the TGU 
approximation given in Eq. fl68|) . However, the separation of the stable and unstable planar 
orbits is still very small. Therefore, the Gutzwiller trace formula dramatically overestimates 
their combined amplitudes, as seen from the short-dashed line. 



5 Summary 

We have derived uniform approximations for the semiclassical description of systems with per- 
turbed SU(2) and SO (3) symmetry. Different from the case of bifurcations, where uniform 
approximations are most effectively derived from the expansion of the classical action in phase- 
space into normal forms we use the group integral representation of the unperturbed trace 
formula || and the classical perturbation theory |14| as starting points for our development. 



In terms of analytical trace formulae we can interpolate smoothly from the integrable limits 
to the limits where their symmetries are broken and the asymptotic Gutzwiller trace formulae 
for the leading periodic orbits are reached. In the two-dimensional Henon-Heiles type poten- 
tials with broken SU(2) symmetry, these orbits are isolated, whereas in the axially-symmetric 
quadrupole-deformed cavities with broken SO (3) symmetry, they retain a one-dimensional de- 
generacy corresponding to the rotation about the symmetry axis. In all cases, the few shortest 
periodic orbits that were considered could account quantitatively for the gross-shell effects 
found in the coarse-grained quantum-mechanical level densities. 

Unlike for the breaking of an orbit family with U(l) symmetry into pairs of stable and 
unstable isolated orbits, for which case Tomsovic et al. fl9| have found the universal uniform 
approximation given in Eq. (|TTJ) , our present results for SU(2) and SO (3) breaking are not 
universal. This is due to more complicated scenarios for the breaking of rational tori that 
arise when more degrees of freedom, or higher-dimensional degeneracies of the orbit families, 
are involved. The two examples of the quartic and cubic Henon-Heiles potentials show how 
the symmetry breaking can happen at different orders in the perturbation theory and lead to 
quite different modulation factors; see Eqs. ( |T9|) and fl45|) which were the starting point of our 
development. However, due to the feasibility of the perturbation theory and the generality of 
the scheme for the uniform approximation used in this study, our method can easily be applied 
also to other potentials and different symmetries of the unperturbed system. 

We acknowledge a critical reading of the manuscript by Martin Sieber, and further helpful 
discussions with Ch. Amann, J. Blaschke, S. Creagh, S. Fedotkin, J. Law, and S. Tomsovic. This 
work has been partially supported by Deutsche Forschungsgemeinschaft (grant No. Br 733/9-1) 
and by Deutscher Akademischer Austauschdienst (DAAD). 
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A Analytical results for linear periodic orbits 



For the test of numerical routines that solve the equations of motion and determine the ac- 
tions of periodic orbits, it might be helpful to compare with analytical results where available. 
The straight-line librating periodic orbits in the Henon-Heiles potentials allow for an analytic 
calculation of their actions or periods. We give here the results and a brief sketch of their 
derivations. 



A.l The Ai and A2 orbits in the quartic Henon-Heiles potential 

The potential 

V(r,6) = ;Vr 2 -^r 4 cos(4#), (77) 

has two pairs of linear librating orbits: the Ai orbits oscillating along the x and y axes (8 = 
and 7r/2), and the orbits A 2 oscillating along the diagonals (9 = 7r/4 and 3tt/4). We scale the 
potential with the factor 1/E* = 4a /uj 4 , so that the equations for their classical turning points 
are 

e = 2x 2 T x 4 (78) 

in terms of the scaled energy e ([[7]) and the scaled radial coordinate x = ryfa/uj. The upper 
and lower sign in ( [78] ) holds for the orbits Ai and A 2 , respectively. The four solutions of Eq. 
(78) with the "— " sign are given by ±x\ and ±£2? where 



x x = V 1 - VI -e, x 2 = V 1 + VI -e, (79) 
and ±Xi are the classical turning points of the Ai orbits. Their action is then 

S Al = 8V2 — [ X \/(e - 2x 2 + x 4 ) dx = 8V2 — r\l(x\-x 2 ){xl-x 2 )dx (80) 

UJ Jo U! JO 

which can be expressed in term of the complete elliptic integrals K(t) = F(ir/2,t) and E(£) = 
E(n/2,t) in terms of the quantity 



2 l ~^ e (81) 



x 2 J l + vr^ 



The result is 



S Al (e) = ^ — yjl + Vl^~e \E(t) - vT^K(t)l . (e < 1) (82) 

At e = 1 we get simply S Al (l) = I6V2 E*/{3u). 

For the A2 orbits, the solutions of Eq. ( [78] ) with the "+" sign are given by ±x± and ±ix 2 , 
where now, 



x x = Vvl + e-l, x 2 = VVl + e + 1, (83) 
and ±Xi are again the classical turning points. The action of the A 2 orbits is then 

S A2 = 8V2 — ( X \l{xl + x 2 ){xl-x 2 )dx (84) 
UJ Jo 
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which becomes 



3 uj 

in terms of the quantity 



SaM = l 4 — (1 + e) 1/4 [(l + VTT~e) K(k) - 2 E(«)l , (85) 



a; 



y/T+e-1 



(86) 



Taylor expansion of Eqs. (^,^) in powers of e leads to the result given in Eq. (|2~5|) . 

A. 2 The A orbit in the standard Henon-Heiles potential 

The potential 

V(r,9) = -wV- -r 3 cos(3#) (87) 
2 3 

has linear librating orbits A oscillating along the symmetry axes (6 = 0, 27r/3 and An/ 3). We 
scale the potential with the factor 1/E* — 6a 2 /uj 6 , so that the equation for the classical turning 
points is 

e = 3x 2 - 2x 3 (88) 

in terms of the scaled energy e = E/E* and the scaled radial coordinate x = ra/uj 2 . The real 
solutions for this cubic equation for e < 1 are, with x\ < x 2 < x 3 , 

xi = - - cos(7r/3 - 0/3) , x 2 = - - cos(7r/3 + 0/3) , x 3 = - + cos(0/3) , (89) 

Zi Zi JLi 

where 

cos0=l-2e. (90) 

The action of the A orbit is then 

Sa_ = aV3 — / V e — 3x 2 + 2x 3 dx = AVQ — / J (x — Xi)(x 2 — x)(x 3 — x) dx , (91) 

but we could not find an analytical expression for this integral. Instead, we calculate the period 
Ta = dSA/dE 

T A = — / dx (92) 

U Jxi J (x — Xi)(x2 — x)(xs — X) 

which again can be expressed in terms of a complete elliptic integral by 

TA(e) = ^-==K(g), ( e < 1) (93) 

where 

5 =(^). (94) 
\x 3 - x\) 

(Note that Ta diverges at e = 1.) Expansion of (|93D in powers of e and integrating over the 
energy E leads to the result given in Eq. (^7|). 
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